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Abstract 

Nuclear scissors modes are considered in the frame of Wigner function moments 
method generalized to take into account spin degrees of freedom and pair correlations 
simultaneously. A new source of nuclear magnetism, connected with counter-rotation of 
spins up and down around the symmetry axis (hidden angular momenta), is discovered. 

Its inclusion into the theory allows one to improve substantially the agreement with exper¬ 
imental data in the description of energies and transition probabilities of scissors modes 
in rare earth nuclei. 
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1 Introduction 

The nuclear scissors mode was predicted [l]-[4] as a counter-rotation of protons against neutrons 
in deformed nuclei. However, its collectivity turned out to be small. From RPA results which 
were in qualitative agreement with experiment, it was even questioned whether this mode is 
collective at all [5, 6]. Purely phenomenological models (such as, e.g., the two rotors model 
[7]) and the sum rule approach [8] did not clear up the situation in this respect. Finally in a 
recent review [9] it is concluded that the scissors mode is ’’weakly collective, but strong on the 
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single-particle scale” and further: ’’The weakly collective scissors mode excitation has become 
an ideal test of models - especially microscopic models - of nuclear vibrations. Most models 
are usually calibrated to reproduce properties of strongly collective excitations (e.g. of J"" = 2+ 
or 3“ states, giant resonances, ...). Weakly-collective phenomena, however, force the models 
to make genuine predictions and the fact that the transitions in question are strong on the 
single-particle scale makes it impossible to dismiss failures as a mere detail, especially in the 
light of the overwhelming experimental evidence for them in many nuclei [10, 11].” 

The Wigner Function Moments (WFM) or phase space moments method turns out to be 
very useful in this situation. On the one hand it is a purely microscopic method, because it is 
based on the Time Dependent Hartree-Fock (TDHF) equation. On the other hand the method 
works with average values (moments) of operators which have a direct relation to the considered 
phenomenon and, thus, make a natural bridge with the macroscopic description. This makes it 
an ideal instrument to describe the basic characteristics (energies and excitation probabilities) 
of collective excitations such as, in particular, the scissors mode. 

Further developments of the WFM method, namely, the switch from TDHF to TDHF- 
Bogoliubov (TDHFB) equations, i.e. taking into account pair correlations, allowed us to im¬ 
prove considerably the quantitative description of the scissors mode [12, 13]: for rare earth 
nuclei the energies were reproduced with ~ 10% accuracy and B(M1) values were reduced by 
about a factor of two with respect to their non superfluid values. However, they remained 
about two times too high with respect to experiment. We have suspected, that the reason of 
this last discrepancy is hidden in the spin degrees of freedom, which were so far ignored by the 
WFM method. 

In a recent paper [14] the WFM method was applied for the first time to solve the TDHF 
equations including spin dynamics. As a first step, only the spin orbit interaction was included 
in the consideration, as the most important one among all possible spin dependent interactions 
because it enters into the mean held. The most remarkable result was the discovery of a new 
type of nuclear collective motion: rotational oscillations of ” spin-up” nucleons with respect of 
’’spin-down” nucleons (the spin scissors mode). It turns out that the experimentally observed 
group of peaks in the energy interval 2-4 MeV corresponds very likely to two different types of 
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motion: the orbital scissors mode and this new kind of mode, i.e. the spin scissors mode. The 
pictorial view of these two intermingled scissors is shown on Fig. 1, which is just the modihcation 
(or generalization) of the classical picture for the orbital scissors (see, for example, [7, 9]). 



Figure 1: Pictorial representation of two intermingled scissors: the orbital (neutrons versus 
protons) scissors + spin (spin-up nucleons versus spin-down nucleons) scissors. Arrows inside 
of ellipses show the direction of spin projections, p - protons, n - neutrons. 

The next step was done in the paper [15], where the influence of the spin-spin interaction 
on the scissors modes was studied. There was hope that, due to spin dependent interactions, 
some part of the force of Ml transitions will be shifted to the energy region of 5-10 MeV (the 
area of a spin-flip resonance), decreasing in such a way the Ml force of scissors. However, 
these expectations were not realised. It turned out that the spin-spin interaction does not 
change the general picture of the positions of excitations described in [14] pushing all levels 
up proportionally to its strength without changing their order. The most interesting result 
concerns the B(M1) values of both scissors - the spin-spin interaction strongly redistributes 
Ml strength in favour of the spin scissors mode practically without changing their summed 
strength. 

In the present work we suggest a generalization of the WFM method which takes into 
account spin degrees of freedom and pair correlations simultaneously. According to our previous 
calculations these two factors, working together, should improve considerably the agreement 
between the theory and experiment in the description of nuclear scissors modes. 
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The paper is organized as follows. In Sec. 2 the TDHFB equations for the 2x2 normal and 
anomalous density matrices are formulated and their Wigner transform is found. In Sec. 3 the 
model Hamiltonian and the mean field are analyzed. In Sec. 4 the collective variables are dehned 
and the respective dynamical equations are derived. In Sec. 5 the choice of parameters and the 
results of calculations of energies and B(M1) values of two scissors modes are discussed. The 
phenomenon of counter-rotating angular momenta with spin up/down, which can be considered 
also as a phenomenon of hidden angular momenta, is analysed in Sec. 6. Results of calculations 
for 26 nuclei in the rare earth region are discussed in Sec. 7. The summary of main results is 
given in the conclusion section. The mathematical details are concentrated in Appendices A, 
B, C, D. 

2 Wigner transformation of TDHFB equations 


The Time-Dependent Hartree-Fock-Bogoliubov (TDHFB) equations in matrix formulation 
are [16, 17] 


with 


ihn = [77, n] 



( 1 ) 

( 2 ) 


The normal density matrix p and Hamiltonian h are hermitian whereas the abnormal density 
k and the pairing gap A are skew symmetric: k) = —k*, = —A*. 

The detailed form of the TDHFB equations is 


ihp = hp — ph — Ak^ -|- kA^, 

—ihp = h*p* — p*h* — A^k -|- k^A, 

—ihk = —hk — kh* -|- A — Ap* — pA, 

—ikn = h*k^ + k^h — A'*^ -|- A'*^p -|- p*A^ (3) 


It is easy to see that the second and fourth equations are complex conjugate to the first and 
third ones respectively. Let us consider their matrix form in coordinate space keeping all spin 
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indices s, s', s": 


ih{r, s|p|r", s") = 

Y, f d^r' ((r,s|/i|r',s')(r',s'|p|r",s") - (r, s|p|r', s')(r', s'|/i|r", s") 

s' 

-(r,s|A|r',s')(r',s'|fi;^|r",s") + (r, s|fi;|r', s') (r, s'| A^", «")) , 
ih{r,s\k\v”,s'') = -(r,s|A|r",s") 

+ Y Jd^^' {jSi s\h\r', s')(r', s'|fi;|r", s") + (r, s|fi:|r', s')(r', s'\h*\r'', s") 

s' 

+ (r,s|A|r',s')(r',s'|p*|r",s") + (r, s|p|r', s')(r', s'|A|r", s")) , 

ih{r, s|p*|r", s") = 

^ J dV (-(r,s|/;,*|r',s')(r',s'|p*|r",s") + (r, s|p*|r', s')(r', s'|/i*|r", s") 

5 ' 

+ {r,s\A^r',s'){r',s'\k\r",s") - (r, s|kV', s') (r', s'| A|r", s")) , 
ih{r,s\K\r'',s") = (r, s| A'*'|r", s") 

k-Y J d^A (—(r, s|/i*|r^ sO(r', s'|fi:'*'|r", s") — (r, s|fi:^|r', s') (r', s'|/i|r", s") 

5 ' 

-(r,s|At|r',s')(r',s'|p|r",s") - (r, s|p*|r', s')(r', s'|A^|r", s")) . 


( 4 ) 


We do not specify the isospin indices in order to make formulae more transparent. They will 
be re-introduced at the end. Let us introduce the more compact notation (r, s|X|r', s') = 

Then the set of TDHFB equations (4) with specihed spin indices reads 


ihpll„ 

ihpH,, 

ihp^„ 

ihkll,, 

ikkil,, 


f ^3 r/i^n tt tt z,tt , iU tt tt z,tt Att I tt AfB ^ 

I U I yitrprpi ^ fj^^f I J ^ 

f -I- 'l 

I tl I fjy.y.1 I iy.fy.1/ ^ ! t y, t fJ y,! y.H fJ yyll t y.t y.t I j ^ 

Jd^r'{h^l,pY, - p^l,hlt„ + h^,p^}'^„ - pti,hti„ - 

-A];j:„ + Jd^r' (hll,tAXn + KXh*X" + ^Xp*X" + pll'Kt") , 

-A;t))„ -h f d^r' {^hXXt" + Xl'd*Vr" + ^tl'P*Vr" + ■ 


( 5 ) 


This set of equations must be complemented by the complex conjugated equations. Writing 
these equations, we neglected the diagonal matrix elements in spin, K^t and A^*/. It is shown in 
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Appendix A that such approximation works very well in the case of monopole pairing considered 
here. 


We will work with the Wigner transform [17] of equations (5). The relevant mathematical 
details can be found in [12], The most essential relations are outlined in Appendix B. /^From 
now on, we will not write out the coordinate dependence (r, p) of all functions in order to make 
the formulae more transparent. The Wigner transform of (5) can be written as 

ihf^ 




- f /"}} + f {{ft }} + -tA- - A^- 

+ f {k. A-} - |{A.k-} - |{{k. A-}} + ^{{A.k-}} + 


ihf^ = 




?{{ft*t /"}} + ?{{/". A''*}} + A-s - s-A 


+ f {A-,k} - |{v, A} - ^{{A-,k}} + ^{{V, A}} + 

/«(;,n _ ,,A) + |{(,,n + yti} _ ;«}} 

,,-tHf-n _ fU) + (;tt + yU)} + (;n _ ;U)}} + 


ihf^^ = 


ihk = 


ihk* = 


/«(/,« _ ;,tt) + y«} _ ;«}} 

A"(/" - !") + |{A*t (/" + /")} + f {{A". (/'^ - f t}} + .... 

H + Z{(A^^ “ AA)^ ,;} — ft}} 

2 8 

A a" + ft + |{(/Af - ft, A} - ^ {{(f f + /A), A}} - A + .... 

-ft"(fc" + A*t + _ /!«), K*} + ^{{(/fAf + A"), K*}} 

2 8 

At/" + f A) + *{(!" - /"), A-} + ^{{(f A + /"), A-}} + A- + 


( 6 ) 


where the functions h, /, A, and n are the Wigner transforms of h, p, A, and k, respectively, 
/(r, p) = /(r, — p), {/, g} is the Poisson bracket of the functions /(r, p) and p(r, p) and 
is their double Poisson bracket; the dots stand for terms proportional to higher powers of h. This 
set of equations must be complemented by the dynamical equations for k, k*. 
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They are obtained by the change p —?> —p in arguments of functions and Poisson brackets. So, 
in reality we deal with the set of twelve equations. We introduced the notation k = and 
A = Symmetry properties of matrices k, A and the properties of their Wigner transforms 
(see Appendix B) allow one to replace the functions fi;'*'^(r,p) and A'*'^(r,p) by the functions 
p) and A^'*'(r,p). 

Following the paper [14] we will write above equations in terms of spin-scalar 

f+ ^ jn + fU 


and spin-vector 


f- = f' - /" 


functions. Furthermore, it is useful to rewrite the obtained equations in terms of even and 
odd functions /g = ^(/ + /) and fo = |(/ — /) and real and imaginary parts of k, and A; 
= 1(k + K*), K* = - K*), A’’ = |(A + A*), A* = T(A - A*). We have 


+ [{fci. /in+{fcit /in+{fct. /in+{/it /in; 

+ 4, ([«■ a:] - [<a:] + |<a;] - [k'a;]) + 

is/o* = y [{y./t} + {y./t} + {ft/./o“} + {ft/.//} 

+ »ft [{fti. /in+(ftit /in+(ftit /in+(ftit /in; 

+ 2,h (K, a:} + {<, A}} + K, a;} + {<, a;}) + 

+ hl^fo^ 

ifl r 1 

+ Y {ht,fe} + {h^,fo} + {h^,f^} + {h^,f^} 

+ 2,n ({k/ a/} + {<, A}} + K, a:} + {<, a:}) 

- Y n"}} + {{ftl*. fi'}} - {{'P. /ini - {{ftf. /in}] + -. 

*// = - /if/" - /if/") 

+ 4i ([«■ Af - [K^Af + [k-A f - |<A}]) 

if) r ~\ 

+ Y {K1 fo} + {ht 1 fe) + {ho 1 ft) + {he 1 ft} 
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J L 


{{hi\ m}+{{hi\ m} - {{hi\ m} - {{hi\ m 


+ 


+ I [{*:*. /„n+ {ii\ fn + {ht, /„"}+{!>:, e) 




ihf^^ = - 
ih 
~2 


= 


= - 
ih 


{{ft?, fe)} + {{ft?, fo)} - {{ftr. /?}} - {{ft„-. /?}} 

Kff + Kfo'' - ft?/7 - ft?/o“ 

+ f [{ft?. /„*} + {ft?. /?} + {ft? /„*? + {ft? /?}[ 

{{ft?. ??} + {{ft?, fo}} - {{ft;. /?}} - {{ft;, fi'}} 
ft;/? + ft;/? - ft?/; - ft?/; 

+ I [{ft?. /?} + {ft?. /„? + {ft? /?} + {ft? /?}[ 

+ f [{{ft?. /;}} + {{ft?. /;}} - {{ft;. /?» - {{ft;. /?»[ 

Kfo^ + Kff - h^efo - h^of. 


+ ... 


+ ... 


+ ... 


+ ^ 1 fe } + fo} + {hf ) fe'^} + {^o ) 


2 


/-}}+{{/^f, /-}} - {{K, m] - {{h-, jf}} 


+ ... 


ihkl 


ihkh 


ihkl 


ihk\ 


ih 


= + <] + + 

1 f) 

+ ^[ff K + ff K] + jiff: K + /e", A^} - /A: + ..., 

1 T) 

= AK i^'o + K i^l] ++ 

17 ) 

+ * 1 // Aj, + /- A'j + -{/„+, a; + /-, a:} - iAj, + ..., 

iT) 

= -AK i^l + K 1^1] ++ 

1 f) 

- ?/?a:+/; a:i + ^{/?,a: + /.-,Ay+ ia: + 

if) 

= -i[ht < + K f^l] ++ 

11 ) 

- ?/?a; + /; Aa + -{/„?A- + /,-,Ay + iA; + 


( 7 ) 


The following notation is introduced here: h^ = ± h'*''*', [ab] = ab— ^{{a, &}}, [ab + cd + 

...] = [flfe] T [cd] -? ...j {n, 6 T c, d T ...}■ = {n, &}■ T {c, dj- T _ 
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These twelve equations will be solved by the method of moments in a small amplitude 
approximation. To this end all functions /(r, p, t) and K,{r, p, t) are divided into equilibrium part 
and deviation (variation): /(r, p, t) = /(r, p)eg + Sf{r, p, t), K,{r, p, t) = K,{r, p)eq + p, t). 
Then equations are linearized neglecting qnadratic terms. 

From general argnments one can expect that the phase of A (and of k, since both are 
linked, according to equation (20)) is mnch more relevant than its magnitnde, since the former 
determines the superflnid velocity. After linearization, the phase of A (and of k) is expressed 
by (5A* (and Sk^), while SA"^ (and Sk^) describes oscillations of the magnitnde of A (and of k). 
Let ns therefore assnme that 

(5K’’(r, p) < (5fi;*(r, p). (8) 

This assnmption was explicitly conhrmed in [18] for the case of snperflnid trapped fermionic 
atoms, where it was shown that SA^ is snppressed with respect to <5A* by one order of A/Ep, 
where Ep denotes the Fermi energy. 

The assnmption (8) allows one to neglect all terms containing the variations and 6A^ in 
the equations (7) after their linearization. In this case the ’’small” variations and (5A” will 
not affect the dynamics of the ’’big” variations Sk^ and 5A* . This means that the dynamical 
eqnations for the ’’big” variations can be considered independently from that of the ’’small” 
variations, and we will hnally deal with a set of only ten eqnations. 


3 Model Hamiltonian 


The microscopic Hamiltonian of the model, harmonic oscillator with spin orbit potential plus 
separable quadrupole-quadrupole and spin-spin residnal interactions is given by 


with 


A 

^ = E 

i=l 


P? 1 2 2 1 e 

— + -mo, r. - „1.S. 


+ Hqq + Hss 


2 f Z N 

^ E (-ir i 

H=-2 [ i j 

1 ( Z N 

n..= Ef-^r{xEEA 

/^=-l I i j 


K 

2 

X 


Z N 

E + E 

Z N 

E + E 


92-,.(ri)9'2^(r,), 
S-„(!)S„(j)i5(rj - rj), 


( 9 ) 

( 10 ) 

( 11 ) 
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where N and Z are the numbers of neutrons and protons and are spin matrices [19]: 

3.1 Mean Field 

Let us analyze the mean held generated by this Hamiltonian. 

3 . 1.1 Spin-orbit Potential 

Written in cyclic coordinates, the spin orbit part of the Hamiltonian reads 


his = -V Y. = -v( 

fi=-i \ —b 


1 h 1 h \ 

I 02 

h 1 h 

V2 ''^2 ' 


where [19] 


l, = -hV2YcYo.ruVa, (13) 

iy,a 

cyclic coordinates r_i,ro,ri are dehned in [19], is a Clebsch-Gordan coefficient and 

k = h(roVi - riVo) = ~^ik + ily), k = h(r_iVi - riV_i) = b, 

Li = h(r_iVo - roV_i) = - ky), 

lx z y) , ly X 2)5 k y x') ‘ (l^) 


Matrix elements of his in coordinate space can be obviously written [14] as 

(ri,si|h«,|r2,S2) = -^V [lokkiSs.^Ss^t - Ss.iSs^i) 

+ v/2/_i(ri)5sit<5s2.^ - V2b(ri)(5si.^552tj - r2)- 


(15) 


Their Wigner transform reads [14]: 


h"Yi^: P) = -^V [^o(r, p)(( 5 ,,t< 5 . 2 t - <^sa 5 s 2 ^) 

+ V2l_i{r, - V^b(r, p)<5si.^<5s2t 


( 16 ) 


where k = -'iV^12u,aCYarvPa- 
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3.1.2 Quadrupole-quadrupole interaction 


The contribution of Hqq to the mean held potential is easily found by replacing one of the g 2 /i 
operators by the average value. We have 

(17) 


Here 


= x-Rt + xK . Zh = Ai;(«) = T /<i(P. r)« 2 ,(r)/"+(r, p, t) (18) 

with /(i(p,r) = {2nh)~^ Jd^p Jd^r and r being the isospin index. 


3.1.3 Spin-spin interaction 


The analogous expression for Hgs is found in a standard way [15] with the following result for 
the Wigner transform of the proton mean held: 




+ Xy + {Sstds't - 

Sp 


(19) 


/ (jj 'P / 

(2 The Wigner transform of the neutron mean held 

is obtained from (19) by the obvious change of indices p ^ n. 


3.2 Pair potential 

The Wigner transform of the pair potential (pairing gap) A(r, p) is related to the Wigner 
transform of the anomalous density by [17] 

A(r,p) = - J - p\)K{r,p), (20) 

where v{p) is a Fourier transform of the two-body interaction. We take for the pairing interac¬ 
tion a simple Gaussian of strength Vq and range Vp [17] 

v{p) = (21) 

with /3 = —|Vo|(rpA/F)^ and a = r'^/Ah^. For the values of the parameters, see section 5.1. 
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4 Equations of motion 

Integrating the set of equations (7) over phase space with the weights 


W = {r 0 p}a„, {r 0 r}A„, {p 0 p}a,„ and 1 


( 22 ) 


one gets dynamic equations for the following collective variables: 


= Jd{P, r){r O p}x^^Sf^^{r, p, t), Tlllit) = Jd{p, r){r O r}x^^6f^^{r, p, t), 

= yd(p,r)(5/;^(r,p,t), 


^xi,it) = / dip, r){r 0 p}A^(5<*(r, p, t), TZl^it) = / d(p, r){r 0 rjx^Snl^r, p, t), 


V 


'v(^) = / dip, r){p O p}xf,5Kl^ir, p, t), 


where ?= +, ti, it, and {r 0 r}xf, = 


( 7 , 1 / 


The required expressions for and are 


(23) 


K = —+ ma;V + 12^(-l)^ZJ+(t){r(g)r}2_/, + h;+(r,t) - 
m ^ p 


being the chemical potential of protons (r = p) or neutrons (r = n), 

h; =-hr//o + l/“(r,t), hl^ = -^pl-i + V}^iY,t), hf = ^r//i + V;'^^(r, t), 
where according to (19) 


.h 






=-3yxn+(r,t), 1/ (r, f) = 3y (r, t) + —(r, t), 

= 3yxnj;'^(r,t) + ^xn);'^(r, t), = 3yxnJ^(r,t) + (24) 


and the neutron potentials are obtained by the obvious change of indices p ^ n. Variations 
of these mean helds read: 


Sh+ = 12j:(-l)'‘iZJ,|(«){r0r}2_„ + «';+(r,«), 




where = kSH^; + hSR^+, SRl^it) = Kl+(t) and 


2fi 


''Xfi ’ 


<^^p^(r,t) =-3yX(5n+(r,t), 6 n^ir,t) = J J^^df^ir,p,t). 


(27rh)" 
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Variations of h , and are obtained in a similar way. Variation of the pair potential is 

(5A(r, p,t) = - J “ P P'> t). (25) 

We are interested in the scissors mode with quantum number K'^ = I"*". Therefore, we only 
need the part of dynamic equations with /i = 1. 

It is convenient to rewrite the dynamical equations in terms of isoscalar and isovector 
variables 

KA,. = n. + Kv. % = + = + 

n,^ = ni^-n% = = (26) 

It also is natural to dehne isovector and isoscalar strength constants ki = — R) and Kq = 

\{k+R) connected by the relation = anQ [20]. Then the equations for the neutron and proton 
systems are transformed into isovector and isoscalar ones. Supposing that all equilibrium 
characteristics of the proton system are equal to that of the neutron system one decouples 
isovector and isoscalar equations. This approximations looks rather crude, nevertheless the 
possible corrections to it are very small, being of the order The integration yields the 

following set of equations for isovector variables: 


£21 = —’P 21 ~ — iy/SanoRoo + V^(l + tt)^o-R2o ~ Ai + 2^22 + ; 

T7Z ^ 



13 





10 


= -n 


pit 

n 


21 




22 


'^20 — 


2v^L 


iCp + hF^^' 


= -V 


£ii + V^/^io 


= -c 


^21 ~ 


m 

2 


21 




2 L 


7^.21 + 27^.22 + V^7^2C) 


— 

> 21 — 


_£- _ ifil']Z+ 

^rll _ ihlTZ+ 

^22 ''"'c^'^211 

m 2 




m 


20 




21 ) 


mo;^ + v^ko-R2o + 6v^fi;o7?20'^n ~ ”^21 + 27^22 + '^'^20 


3^3 




*4-1 *4-2 


(A- 


T^in 


22 


pp 


20 


/’it 

-^00 


-^l^olC(r')P2i, 


^2"l = 


mu^ + V/6K0-R20 Ai + Sv^Ko-R^O'^ii ~ 6 v^koLio(®ii)"^Ji ~ i^^ 7 T^ 2 i 


3^3 




*4i*42 


(4-1 — 4 . 2 ) F 21 + (4i + 42 ) F 
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h 


2m u'^ - aVGkqRII - ^^h^x . 

I A.2 


cli - ih^v. 


21 ) 


2 

./3, 


2m+ 4v/6Koi?2o -^^o + Sv/S/toi^so'C^ - -Vp 


a/3 2 -^2 

' 2 ^ ^4i42 


(4i — 242 ) (4i + 42 ) Z^oo 


= -vil - m + Ad^KoRll n^o 


m 


2^3 




X -fjh- Ixh 


J [V^2 


^ ^ I ^ 

42 4 


)tt 

-20 


T^tt _ 
'^00 — 


pit 

' 00 


7^21 

’P 21 

di 

Rii 


m 


rll 

-^00) 




-2mdC\l + sVSKoi^ao -^20 + 


_L '\ /’-^t 

42 4i 


^ (^yaKo^4 + Ao(r') - 


n 


21 ) 


■-Ao(r')7^2i + QhaKoFoTZ 


+ 

21 ) 


Ao(r')A 


21 ) 




00 


^ ± ] £^t 

42 4i ' 


20 


(27) 
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where A = i ?20 - ^oo = ^ (l + |<5) , A = i? 2 o/v^ + ^oo = (l “ 1 ^), «-i = 

^ fl-i2/3)6Y^^ , „ ^l-(2/3)5V'^' • ^ n- •. u u-u.u 

ai = Kq \ - ,, , ■ , and an = Kn \ - ,, , . ^ are semiaxes oi ellipsoid by which the 

'’\l + {A/3)6) ^\l + {A/3)6) 

shape of nucleus is approximated, 5 - deformation parameter, Rq = 1.274^/^ fm. 


/i = ^ y drr‘^ (~g~^ ’ J drr‘^n'^{ry, A(r) = 71+ + = 

0 V / n 


1 + e" 


/Co = / (i(r, p)7Co(r, p),/C 4 = / (/(r, p)r^7Co(r, p). The functions 7Co(^"0) 

Ip^ir') are discussed in the next section and are demonstrated in Appendix D. Deriving these 
equations we neglected double Poisson brackets containing k or A, which are the quantum 
corrections to pair correlations. The isoscalar set of equations is easily obtained from (27) by 
taking a = 1, replacing x —x and putting the marks ’’bar” above all variables. 


5 Results of calculations 

The set of equations (27) coincides with the set of equations (27) of the paper [15] in the limit 
of zero pairing, i.e. if to omit the last four equations and to neglect the contributions from 
pairing in the dynamical equations for the variables A^i) V 21 ■ On the other hand, 

the dynamical equations for 7^21 and 77.21 and the contribution from pairing in the dynamical 
equation for V 21 are exactly the same as the ones in the paper [13]. Only the dynamical 
equations for £ 21 , and the contributions from pairing in dynamical equations for £ 21 ; 
are completely new. 

Imposing the time evolution via for all variables one transforms (27) into a set of 

algebraic equations. It contains 23 equations. To hnd the eigenvalues we construct the 23x23 
determinant and seek (numerically) for its zeros. We hnd seven roots with exactly E=0 and 
16 roots which are non zero: eight positive ones and eight negative ones (situation is exactly 
same as with RPA; see [21] for connection of WFM and RPA). In this paper we consider only 
the two lowest roots corresponding to the orbital and spin scissors. The qualitative picture of 
high lying modes remains practically without any changes in comparison with [15]. 

Seven integrals of motion corresponding to Goldstone modes (zero roots) can be found 
analytically. They are written out in the Appendix C. The interpretation of some of them has 
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been found in [15], whereas the interpretation of the remaining ones seems not to be obvious. 


5.1 Choice of parameters 

• Following our previous publications [20, 21] we take for the isoscalar strength constant of the 
quadrupole-quadrupole residual interaction kq the self consistent value [ 22 ] kq = with 

Qoo = R = roAR^, ro = 1.2 fm, - ISyR], huo = MeV. 

• The equations (27) contain the functions Ao(r') = Ae5(r',pi?(r')), 

Ip^{r') = Ip^{r',pF{r')) and ko(^ 0 = depending on the radius r' and the local Fermi 

momentum prir') (see Fig. 2 ). The value of r' is not hxed by the theory and can be used 



Figure 2: The pair held (gap) Ao(r), the function A = \Vo\Ip^{r) and the nuclear density n(r) 
as the functions of radius r. The solid lines - calculations without the spin-spin interaction 1/ss, 
the dashed lines - 17,^ is included. 

as the htting parameter. We have found in our previous paper [13] that the best agreement 
of calculated results with experimental data is achieved at the point r' where the function 
Ip^{r',pF{r')) has its maximum. Nevertheless, to get rid off the htting parameter, we use the 
averaged values of these functions: Aq = /(irno(r)Ao(r,pp’(r))/A, etc. The gap A(r,pj?(r)), 
as well as the integrals Ip^{r,pF{r)), /C 4 and JCq, were calculated with the help of the semiclas- 
sical formulae for K{r, p) and A(r, p) (see Appendix D), a Gaussian being used for the pairing 
interaction with = 1.9 fm and Vq = 25 MeV [17]. Those values reproduce usual nuclear 
pairing gaps. 
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• The used spin-spin interaction is repulsive, the values of its strength constants being taken 
from the paper [23], where the notation x = Kg/A, X = QX was introduced. The constants 
were extracted by the authors of [23] from Skyrme forces following the standard procedure, 
the residual interaction being dehned in terms of second derivatives of the Hamiltonian density 
H{p) with respect to the one-body densities p. Different variants of Skyrme forces produce 
different strength constants of spin-spin interaction. The most consistent results are obtained 
with SGI, SG2 [24] and Sk3 [25] forces. To compare theoretical results with experiment the 
authors of [23] preferred to use the force SG2. Nevertheless they have noticed that ”As is 
well known, the energy splitting of the HF states around the Fermi level is too large. This 
has an effect on the spin Ml distributions that can be roughly compensated by reducing the 
Kg value”. According to this remark they changed the original self-consistent SG2 parameters 
from Kg = 88 MeV, q = —0.95 to Kg = 50 MeV, q = —1. It was found that this modihed 
set of parameters gives better agreement with experiment for some nuclei in the description of 
spin-flip resonance. So we will use Kg = 50 MeV and q = —1. 

• Our calculations without pairing [15] have shown that the results are strongly dependent 
on the values of the strength constants of the spin-spin interaction. The natural question 
arises: how sensitive are they to the strength of the spin-orbital potential? The results of the 
demonstrative calculations are shown in Fig. 3. 

The Ml strengths were computed using effective spin giromagnetic factors = 0.7 
One observes a rather strong dependence of the results on the value of p: the splitting AE and 
the Ml strength of the spin scissors grow with increasing p, the B{M1) of the orbital scissors 
being decreased. At some critical point pc the Ml strength of the spin scissors becomes bigger 
than that of the orbital scissors. The inclusion of the spin-spin interaction does not change the 
qualitative picture, as well as the inclusion of pair correlations (see Fig. 3). 

What value of p to use? Accidentally, the choice of p in our previous papers [14, 15] was 
not very realistic. The main purpose of the hrst paper was the introduction of spin degrees 
of freedom into the WFM method, and the aim of the second paper was to study the influ¬ 
ence of spin-spin forces on both scissors - we did not worry much about the comparison with 
experiment. Now, both preliminary aims being achieved, one can think about the agreement 
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Figure 3: The energies E and i?(Ml)-factors as a functions of the spin-orbital strength constant 
T]. Left panel; solid lines - without the spin-spin interaction Vgs, dashed lines - Vgs is included. 
Right panel: The same as in left panel with pair correlations included. 

with experimental data, therefore the precise choice of the model parameters becomes impor¬ 
tant. Of course, we could try to choose rj according to the standard requirement of the best 
agreement with experiment. However, in reality we are not absolutely free in our choice. It 
turns out that we are already restricted by the other constraints. As a matter of fact we work 
with the Nilsson potential, parameters of which are very well known. Really, the mean held 
of our model (9) is the deformed harmonic oscillator with the spin-orbit potential, the Nilsson 
term being neglected because it generates the fourth order moments and, anyway, they are 
probably not of great importance. In the original paper [26] Nilsson took the spin-orbit strength 
constant k^us = 0.05 for rare earth nuclei. Later the best value of k^hs for rare earth nuclei 
was established [17] to be 0.0637. For actinides there were established different values of knus 
for neutrons (0.0635) and protons (0.0577). The numbers knus = 0.0637, knus = 0.05 and 
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i^Niis = 0.024 (corresponding to rj = 0.36 used in our previous calculations [14, 15]) are marked 
on Figs 3, 5 by the dotted vertical lines. Of course we will use the conventional [17] parameters 
of the Nilsson potential and from now on we will speak only about the Nilsson [26] spin-orbital 
strength parameter knus, which is connected with t] by the relation rj = 2huK,Niis- 


5.2 Discussion and interpretation of results 

The energies and excitation probabilities of orbital and spin scissors modes obtained by the 
solution of the isovector set of equations (27) are displayed in the Table 1. 


Table 1: Scissors modes energies Esc and transition probabilities B{Ml)sc- 


i64gr 

Esc, MeV 

B{Ml)sc, 



Ks = 0 

Ks = 50 

Ks = 92 

Ks = 0 

Ks = 50 

Ks = 92 

spin 

Ao = 0 

1.40 

1.60 

1.73 

5.38 

6.23 

6.79 

scissors 

Ao 7^ 0 

2.72 

2.75 

2.77 

3.93 

4.79 

5.44 

orbital 

Ao = 0 

2.57 

2.69 

2.78 

1.74 

1.59 

1.50 

scissors 

Ao 7^ 0 

3.49 

3.51 

3.52 

1.74 

1.51 

1.35 


There are results of calculations with three values of the spin-spin strength constant and 
two values of Aq. As it was expected the energies of both scissors increased approximately 
by 1 Mev after inclusion of pairing. The behaviour of transition probabilities turned out less 
predictable. The B{M1) value of the spin scissors decreased approximately by 1.5 whereas 
B{M1) value of the orbital scissors turned out practically insensitive to the inclusion of pair 
correlations. 

We can compare the summed B{Ml)-s = B{Ml)or + B{Ml)sp values and the centroid of 
both scissors energies 

Been = [EorB{Ml)or + EspB{Ml)sp\/B{Ml)j: 

with the results of the paper [13] where no spin degrees of freedom had been considered and 
with the experimental data. The respective results are shown in the Table 2. 

It is seen that the inclusion of spin degrees of freedom in the WFM method does not change 
markedly our results (in comparison with previous ones [13]). Of course, the energy changed in 
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Table 2: Scissors modes energy centroid Ecen and summarized transition probabilities B{M1)y,- 


The experimental values of -Ecen, and are from [27] and references therein. 


i64gr 

Been, MeV 

H(Ml)s, lil 


iLs = 50 

iLs = 92 

[13] 

exp 

Ks = 0 

Ks = 50 

iLs = 92 

[13] 

exp 

Ao = 0 

Ao ^ 0 

1.69 

2.96 

1.82 

2.93 

1.92 

2.92 

2.10 

3.37 

2.90 

7.13 

5.67 

7.82 

6.30 

8.29 

6.79 

9.26 

5.62 

1.45 


the desired direction and now practically coincides with the experimental value (especially in 
the case with spin-spin forces.) However, the situation with the B{M1) values did not change 
(and even become worse in the case with spin-spin forces). Our hope, that spin degrees of 
freedom can improve the situation with the B{M1) values, did not become true: the theory so 
far gives four times bigger values of B{M1) than the experimental ones, exactly as it was the 
case in the paper [13]. 

The result look discouraging. However, a phenomenon, which was missed in our previous 
papers and described in the next section will save the situation. 

6 Counter-rotating angular momenta of spins up/down 
(hidden angular momenta) 

The equilibrium (ground state) orbital angular momentum of any nucleus is composed of two 
equal parts: half of nucleons (protons -|- neutrons) having spin projection up and other half 
having spin projection down. It is known that the huge majority of nuclei have zero angular 
momentum in the ground state. We will show below that as a rule this zero is just the sum 
of two rather big counter directed angular momenta (hidden angular momenta, because they 
are not manifest in the ground state) of the above mentioned two parts of any nucleus. Being 
connected with the spins of nucleons this phenomenon naturally has great influence on all 
nuclear properties connected with the spin, in particular, the spin scissors mode. 

Let us analyze the procedure of linearization of the equations of motion for collective vari¬ 
ables (23). We consider small deviations of the system from equilibrium, so all variables are 
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written as a sum of their equilibrium value plus a small deviation: 


L(t) = L{eq) + C(t), et al. 


Neglecting quadratic deviations one obtains the set of linearized equations for deviations de¬ 
pending on the equilibrium values R^^i^eq) and L^^(eg), which are the input data of the problem. 
In the paper [15] we made the following choice: 


Rt±i{.eq) = i?^±2(eg) = 0, Rto{eq) ^ 0, i?oo(eg) ^ 0, (28) 

Mi(eg) = ^v(eg) = 0, (29) 

(30) 


At hrst glance, this choice looks quite natural. Really, relations (28) follow from the axial 
symmetry of nucleus. Relations (29) are justified by the fact that these quantities should 
be diagonal in spin at equilibrium. The variables T)(^(t) contain the momentum p in their 
dehnition which incited us to suppose zero equilibrium values as well (we will show below that 
it is not true for Lbo because of quantum effects connected with spin). 

The relation R^^{eq) = 0 follows from the shell model considerations: the nucleons with 
spin projection ”up” and ’’down” are sitting in pairs on the same levels, therefore all average 
properties of the ’’spin up” part of nucleus must be identical to that of the ’’spin down” part. 
However, the careful analysis shows that being undoubtedly true for variables this 

TT I I 

statement turns out erroneous for variables L[q, L|q. Let us demonstrate it. By dehnition 

Lxtit) = J ^ {r,p,t) = J 0 R'}x^^, (31) 

where 


‘) = p. i) = ■ I)'’*'' + (32) 

is the i-th component of the nuclear current. In the last relation the dehnition [17] of Wigner 
function is used. Performing the integration over p one hnds: 

t) = ihjd^q[-^S{q)]p{r + |, s; r - s'; t) 


(Rp 


= -ih Jd^q6{q)-^p{r + ^,s-,r - |,s';t) 
%Ti 

= -y [(Vli - V 2 i)p(ri, s; r 2 , s'; t)]ri=r2=r, 


(33) 
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where ri = r + |, r 2 = r — 2. The density matrix of the ground state nucleus is dehned [17] as 

p(ri, s; r 2 , s'; t) = ^ (ras'), (34) 

V 

where vl are occupation numbers and (pi, are single particle wave functions. For the sake of 
simplicity we will consider the case of spherical symmetry. Then u = nljm and 

0niim(r, s) = Unijir) ^ (p)xiais), (35) 

A,a 


ih 


jfi^) = ■0;:(r,s') -0^(r,s) ■ Vi0;:(r,s'); 


ih 


= -y i; H v.y,;,|x■ ,(i>)x■,,(/)■ 


(36) 

(37) 


nljm A,a,A',a' 

Inserting this expression into (31) one hnds: 

ih 


ro(eg) = -y E 


A(T,A'(t' 

® v}i„y,A - yA{r ® v},„yj,] 


nljm 

'o2 


i: c'iriAC'i 7 iA.xy(i>)xy.(s') dAiziiiiyAWA-YMi,] 

^ nljm Aa,A'a' ’2 ’2 J 


^T.<ijm E c'/;;i^c'/™i^,xi,(s)xi,,(s')(a + a')5a,a' 

nljm 


Aa,A'a' 

i 


= 7 x E <ijmT.^ XlAs)XlaW)- (38) 

Here the dehnition = —h^/2{r ® V}i^, formula loYi\ = ATJa and normalization of functions 
T^nij were used. Remembering the dehnition of the spin function Xio-(s) = ^o-,* we get hnally: 

i 

7h 

nljm A ^ ^ y ^ nljm 

Now, with the help of analytic expressions for Clebsh-Gordan coefficients one obtains the hnal 
expressions 


LToieq) = ^J2<ijmJ2^[Ciris) <ijm [C^Z_sls) ( 39 ) 


Lioieq) = ^E 




E 

nl \rn=—^l+^^ 


2 / + ^ + ?n 


^nljd'm 


21 + 1 


E 

m=-(l-l] 


2 I + I — m 

Vnlj-m 2/ + 1 


m 


, ( 40 ) 
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L\i{eq) = 


72 


nl 


i+h 


E 

Lm=-(/+i) 


/ + I — m 
‘ 2 / + 1 




+ E 

m=— (I— 


„,2 I + I + rn 

2/ + 1 


m+-), (41) 


where the notation E = / ± 4 is introduced. Replacing in (40) m by —m we find that 


4o(eg) = -L^oieq). 


11/ 


(42) 


By definition (23) L^Q{eq) = L^iQ{eq) ± L\^{eq). Combining linearly (40) and (41) one finds: 


rll/ 


Ltoieq) = ^E 

V2 nl 


V 

/ . ^nlj+m 

Lm=-(«+i) 


2 / 


'-5 


2 / + 1 


m 


+ E 


,,2 _ 

^ ^b--2/ + l 


2 /+ 2 




m 


(43) 


v2 


'+1 

E 


1^, 




.™=-('+l) 


2m2 - / - 4 
2 / + 1 


E 

m=— (I— 




2m‘^ + 1 + 1 
2 / + 1 


(44) 


These formulas are valid for spherical nuclei. However, with the scissors and spin-scissors 
modes, we are considering deformed nuclei. For the sake of the discussion, let us consider the 
case of infinitesimally small deformation, when one can continue to use formulae (43, 44). Now 
only levels with quantum numbers +m are degenerate. According to, for example, the Nilsson 
scheme [26] nucleons will occupy pairwise precisely those levels which leads to the zero value 
of Lfoieq). 

What about LiQ{eq)7 It only enters (27) in the equation for V^i- Let us analyze the 
structure of formula (44) considering for the sake of simplicity the case without pairing. Two 
sums over m (let us note them Si and S 2 ) represent two spin-orbital partners: in the first 
sum the summation goes over levels of the lower partner {j = I + ^) and in the second sum 
- over levels of the higher partner (j = / — 4). The values of both sums depend naturally on 
the values of occupation numbers Unijm = 0,1. There are three possibilities. The first one is 
trivial: if all levels of both spin-orbital partners are disposed above the Fermi surface, then 
the respective occupation numbers Unijm = 0 and both sums are equal to zero identically. The 
second possibility: all levels of both spin-orbital partners are disposed below the Fermi surface. 
Then all respective occupation numbers Unij+m = nnij-m = 1- The elementary analytical 
calculation (for arbitrary /) shows that in this case the two sums in (44) exactly compensate 
each other, i.e. Si -|- S 2 = 0. The most interesting is the third possibility, when one part 
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of levels of two spin-orbital partners is disposed below the Fermi surface and another part is 
disposed above it. In this case the compensation does not happen and one gets Si -|- E 2 7 ^ 0 
what leads to Lho(eg) 7 ^ 0. In the case of pairing, things are not so sharply separated and 
Lio{eq) has always a hnite value. However, the modihcations with respect to mean held are 
very small. 

Let us illustrate the above analysis by the example of ^®^Er (protons). Its deformation is 
6 = 0.25 (e = 0.26) and Z= 68 . Looking on the Nilsson scheme (for example. Fig.1.5 of [16] or 
Fig. 2.21c of [17]) one easily hnds, that only three pairs of spin-orbital partners give a nonzero 
contribution to L 7 o(eg). They are: N = 4, ^ 5/2 — <^ 3/2 (two levels of ^ 5/2 are below the Fermi 
surface, all the rest - above); N = 4:,gg/2 — 97/2 (one level of gY/2 is above the Fermi surface, 
all the rest - below); N = 5, hii /2 — hg /2 (four levels of hii /2 are below the Fermi surface, all 
the rest - above). It is possible to make the crude evaluation of LiQ^eq) using the quantum 
numbers indicated on Fig.1.5 of [16] or Fig. 2.21c of [17]). The result turns out rather close 
to the exact one, computed with the help of formulas (31,36) and Nilsson wave functions. The 
inhuence of pair correlations is very small. 


(a) Equilibrium 



(b) Spin-scissors 

U TT 



Figure 4: (a) Protons with spins t (up) and I (down) having nonzero orbital angular momenta 
at equilibrium, (b) Protons from Fig. (a) vibrating against one-another. 

Indeed, from the dehnitions (31) and (38) one can see that Lll{eq) is just the average value 
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of the z-component of the orbital angular momentum of nucleons with the spin projection 
s or —So, the ground state nucleus consists of two equal parts having nonzero angular 
momenta with opposite directions, which compensate each other resulting in the zero total 
angular momentum. This is graphically depicted in Fig. 4(a). 

On the other hand, when the opposite angular momenta become tilted, one excites the 
system and the opposite angular momenta are vibrating with a tilting angle, see Fig. 4(b). 
Actually the two opposite angular momenta are oscillating, one in the opposite sense of the 
other. It is rather obvious from Fig. 1 that these tilted vibrations happen separately in each 
of the neutron and proton lobes. These spin-up against spin-down motions certainly influence 
the excitation of the spin scissors mode. So, classically speaking the proton and neutron 
parts of the ground state nucleus consist each of two identical gyroscopes rotating in opposite 
directions. One knows that it is very difficult to deviate gyroscope from an equilibrium. So one 
can expect, that the probability to force two gyroscopes to oscillate as scissors (spin scissors) 
should be small. This picture is conhrmed in the next section. 

7 Results of calculations continued 

We made the calculations taking into account the non zero value of L^Q^eq) (which was computed 
according to formulas (31,36) and Nilsson wave functions). The results are shown on Fig. 5. 

They demonstrate (in comparison with Fig. 3) the strong influence of the spin-up vs spin- 
down angular momenta on the spin scissors mode, whose B(M1) value is strongly decreasing 
with increasing Kp^ns- The B(M1) value of the orbital scissors also is reduced, but not so 
much, the value of the reduction being practically independent on unus- The influence of 
Lbo(eg) on the energies of both scissors is negligible, leading to the small increase of their 
splitting. Now the energy centroid of both scissors and their summed B(M1) value at kmus = 
0.0637 are Been = 2.97 MeV and B{M1)y, = 3.7/i^. The general agreement with experiment 
becomes considerably better (compare with Table 2), though the theoretical value of i?(Ml)s 
still exceeds the experimental one approximately 2.5 times. However, as we will see, the case of 
i64Er niay imply a quite particular situation (or even a problem with the experimental value). 

The results of systematic calculations for rare-earth nuclei are presented in Tables 3 and 4 
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Figure 5: The energies E and i?(Ml)-factors as a functions of the spin-orbital strength constant 
unus- The dashed lines - calculations without Lho(eg), the solid lines - Li^^eq) are taken into 
account. 14 ^ and pairing are included. 

and desplayed in Fig. 6. Table 3 contains the results for well deformed nuclei with 6 > 0.18. It 
is easy to see that the overall (general) agreement of theoretical results with experimental data 
is substantially improved (in comparison with our previous calculations [13]). 

The results of calculations for two groups (’Tight” and ’’heavy”) of weakly deformed nuclei 
with deformations 0.14 < 6 < 0.17 are shown in the Table 4. They require some discussion, 
because of the self-consistency problem. These two groups of nuclei are transitional between well 
deformed and spherical nuclei. Systematic calculations of equilibrium deformations [16] predict 

= 0.0 for ±0.1 for ^^®Nd, 0.15 or -0.12 for ^®°Sm, 0.1 or -0.14 for ^®°Os and -0.1 for 

^®^Os, whereas their experimental values are 6eq = 0.14, 0.17, 0.16, 0.15 and 0.14 respectively. 
As one sees, the discrepancy between theoretical and experimental 6eq is large. Uncertain signs 
of theoretical equilibrium deformations are connected with very small (~0.1-0.2 MeV) difference 
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Table 3: Scissors modes energy centroid E^en and summarized transition probabilities B{M1)y,. 
Parameters: knus = 0.0637, Vq = 25 {Vq = 27 for I82,i84,i86yy^_ experimental values of Been, 
S, and are from [27] and references therein. 


Nuclei 

5 


-^ceri: 

MeV 



B{M1) 

/,2 

S, 


exp 

present 

[13] 

A = 0 

exp 

present 

[13] 

A = 0 

150Nd 

0.22 

3.04 

2.88 

3.44 

1.92 

1.61 

1.64 

4.17 

7.26 

I52gm 

0.24 

2.99 

2.99 

3.46 

2.02 

2.26 

2.50 

4.68 

7.81 

i54gm 

0.26 

3.20 

3.10 

3.57 

2.17 

2.18 

3.34 

5.42 

8.65 

156Gd 

0.26 

3.06 

3.09 

3.60 

2.16 

2.73 

3.44 

5.42 

8.76 

158Gd 

0.26 

3.14 

3.09 

3.60 

2.19 

3.39 

3.52 

5.72 

9.12 

leoGd 

0.27 

3.18 

3.14 

3.61 

2.21 

2.97 

4.02 

5.90 

9.38 

leopy 

0.26 

2.87 

3.08 

3.59 

2.13 

2.42 

3.60 

5.53 

9.03 


0.26 

2.96 

3.07 

3.61 

2.14 

2.49 

3.69 

5.66 

9.25 

le^Dy 

0.26 

3.14 

3.07 

3.60 

2.17 

3.18 

3.78 

5.95 

9.59 

i64Er 

0.25 

2.90 

3.01 

3.57 

2.10 

1.45 

3.39 

5.62 

9.26 

166 Er 

0.26 

2.96 

3.06 

3.53 

2.13 

2.67 

3.86 

5.96 

9.59 

168 Er 

0.26 

3.21 

3.06 

3.53 

2.10 

2.82 

3.95 

5.95 

9.67 

i70Er 

0.26 

3.22 

3.05 

3.57 

2.09 

2.63 

4.03 

5.91 

9.79 

i^^Yb 

0.25 

3.03 

2.99 

3.55 

2.05 

1.94 

3.72 

5.84 

9.79 

i^^Yb 

0.25 

3.15 

2.98 

3.47 

2.02 

2.70 

3.80 

5.89 

9.82 

176Yb 

0.24 

2.96 

2.92 

3.45 

1.94 

2.66 

3.46 

5.54 

9.58 

178Hf 

0.22 

3.11 

2.81 

3.43 

1.79 

2.04 

2.67 

4.86 

9.00 

180Hf 

0.22 

2.95 

2.81 

3.36 

1.76 

1.61 

2.69 

4.85 

8.97 

182yy 

0.20 

3.10 

3.28 

3.30 

1.63 

1.65 

2.05 

4.31 

8.43 


0.19 

3.31 

3.24 

3.28 

1.55 

1.12 

1.72 

3.97 

8.14 

186yY 

0.18 

3.20 

3.19 

3.26 

1.49 

0.82 

1.40 

3.76 

7.95 
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Table 4: Scissors modes energy centroid Ecen and summarized transition probabilities 


Parameters: k = 0.05 {k = 0.0637 for I82,i84,i86yy^^ 27 


Nuclei 

5 


E 

-^cen? 

MeV 



B{M1) 

/,2 

S, h-AT 


exp 

present 

[13] 

A = 0 

exp 

present 

[13] 

A = 0 

i34Ba 

0.14 

2.99 

3.04 

3.09 

1.28 

0.56 

0.68 

1.67 

3.90 

148Nd 

0.17 

3.37 

3.22 

3.18 

1.48 

0.78 

1.28 

2.58 

5.39 

i50gm 

0.16 

3.13 

3.17 

3.13 

1.42 

0.92 

1.12 

2.45 

5.26 

182y^ 

0.20 

3.10 

3.28 

3.30 

1.63 

1.65 

2.05 

4.31 

8.43 

I84y^ 

0.19 

3.31 

3.24 

3.28 

1.55 

1.12 

1.72 

3.97 

8.14 

186’^ 

0.18 

3.20 

3.19 

3.26 

1.49 

0.82 

1.40 

3.76 

7.95 

igoQs 

0.15 

2.90 

3.14 

3.12 

1.21 

0.98 

1.38 

2.67 

6.64 

1920s 

0.14 

3.01 

3.11 

3.12 

1.15 

1.04 

1.00 

2.42 

6.37 



A 

Figure 6: The energies Egc and i?(Ml)sc-factors as a function of the mass number A for nuclei 
listed in the Table 3. 
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between the values of deformation energies Sdef at positive and negative 6eq- Even more so, 
the values of deformation energies of these nuclei are very small: Sdef = 0.20, 0.50, 0.80 and 
0.70 MeV for ^^®Nd, ^^°Sm, ^®°Os and ^®^Os respectively. This means that these nuclei are very 
” soft” with respect of /3- or 7 -vibrations and probably they have more complicated equilibrium 
shapes, for example, hexadecapole or octupole deformations in addition to the quadrupole one. 
This means that for the correct description of their dynamical and equilibrium properties it is 
necessary to include higher order Wigner function moments (at least fourth order) in addition 
to the second order ones. In this case it would be natural also to use more complicate mean 
field potentials (for example, the Woods-Saxon one or the potential extracted from some of 
the numerous variants of Skyrme forces) instead of the too simple Nilsson potential. Naturally, 
this will be the subject of further investigations. However, to be sure that the situation with 
these nuclei is not absolutely hopeless, one can try to imitate the properties of the more perfect 
potential by htting parameters of the Nilsson potential. As a matter of fact this potential has 
the single but essential parameter - the spin-orbital strength umis- If turns out that changing 
its value from 0.0637 to 0.05 (the value used by Nilsson in his original paper [26]) is enough 
to obtain the reasonable description of B(M1) factors (see Table 4). To obtain the reasonable 
description of the scissors energies we use the ’’freedom” of choosing the value of the pairing 
interaction constant Vq in (21). It turns out that changing its value from 25 MeV to 27 MeV is 
enough to obtain the satisfactory agreement between the theoretical and experimental values 
of Esc (Table 4). 

The isotopes turn out intermediate between weakly deformed and well deformed 

nuclei: reasonable results are obtained with knus = 0.0637 (as for well deformed) and Vq = 27 
MeV (as for weakly deformed). That is why they appear in both Tables. 

Returning to the group of well deformed nuclei with S > 0.18 (Table 3) it is necessary 
to emphasize that all presented results for these nuclei were obtained without any htting. In 
spite of it the agreement between the theory and experiment can be called excellent for all 
nuclei of this group except two: ^®^Er and ^^^Yb, where the theory overestimates B(M1) values 
approximately two times. However, these two nuclei fall out of the systematics and one can 
suspect, that there the experimental B(M1) values are underestimated. Therefore one can hope. 
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that new experiments will correct the situation with these nuclei, as it happened, for example, 
with 232Th [28], 


8 Conclusion 

The method of Wigner function moments is generalized to take into account spin degrees of 
freedom and pair correlations simultaneously. The inclusion of the spin into the theory allows 
one to discover several new phenomena. One of them, the nuclear spin scissors, was described 
and studied in [14, 15], where some indications on the experimental conhrmation of its existence 
in actinides nuclei are discussed. Another phenomenon, the opposite rotation of spin up/down 
nucleons, or in other words, the phenomenon of hidden angular momenta, is described in this 
paper. Being determined by the spin degrees of freedom this phenomenon has great influence 
on the excitation probability of the spin scissors mode. On the other hand the spin scissors 
B(M1) values and the energies of both, spin and orbital, scissors are very sensitive to the 
action of pair correlations. As a result, these two factors, the spin up/down counter-rotation 
and pairing, working together, improve substantially the agreement between the theory and 
experiment in the description of the energy centroid of two nuclear scissors and their summed 
excitation probability. More precisely, for the hrst time an excellent agreement is achieved for 
well deformed nuclei of the rare earth region with standard values of all possible parameters. 

An excellent agreement is also achieved for weakly deformed (transitional) nuclei of the 
same region by a very modest re-£t of the spin-orbit strength. We suppose that fourth order 
moments and more realistic interactions are required for the adequate description of transitional 
nuclei. However this shall be the subject of future work. 
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According to formula (D.47) of [17] the abnormal density in coordinate representation 
K{r, s] r', s') is connected with the abnormal density in the representation of the harmonic 
oscillator quantum numbers = (<h|a,^a,^/|<h) by the relation 

s; r', s') = (<h|a(r, s)a(r', s')|<l>) = s)^/’^/(r', s')(<l>|a,,a^/|<l>), (A.l) 

where 

u = k,<;, k = n,l,j,\m\, ^ = sign(m) = ±, A;, + = z/, k,— = v, 

^/>p(r,s) = T^7^(r,s). 

T - time reversal operator dehned by formula (XV.85) of [29] 

T KTyK^^ 

where ay is the Pauli matrix and Kq is the complex-conjugation operator. 

According to formula (7.12) of [17] 

Q;,,|<h) = 0 , 

(<h|a,,a,,/|<l>) = ■ (A.2) 

This result means that in accordance with the theorem of Bloch and Messiah we have found 
the basis \v) in which the abnormal density has the canonical form. Therefore the spin 
structure of K^y is 

i^vy = 

or and 0* 

With the help of (A.2) formula (A.l) can be transformed into 

K(r,s;r',s') = s)^k-<;{r , s') 

/c,<7 

= ^M,,n,,[^/’,,(r,s)^/’p(r',s')-^/’p(r,s)^/’,,(r',s')]. (A.4) 

u>0 

that reproduces formula (D.48) of [17]. 


0 

~'^kVk 


0 


(A.3) 
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What is the spin structure of K{r, s-,r', s')7 
Let us consider the spherical case: 


^^(r, s) = Unijir) 2 ^ 1^1 

A,a 






where s) = a ^za(^, function xio-('S) = ^ 0,3 and angular variables 

5 t/V, 2^2 2 

are denoted by f 2 . 

Time reversal: 

rviA = y,i = 

rx..=xi-i, rx._.=-xy ^ rxi, = (-ir^hi-,, 

TEclZ.Mxi. = ECiT:i,«-AX|-A(-i)""-’ = ECiTLi-,«AXiA(-i)- 

A,a ^ A rr ^ A A. 2 


\ —A—cr —; 


^ A,fT ^ A,a ^ 

=Ecy5,iAxiA(-i)'+’-'’-'-'’-^ = EcqpiAX.A(-iy-'*’” 

A,(7 A,cr 


As a result 




that coincides with formula (2.45) of [17]. Formula (A.4) can be rewritten now as 


K(ri,Si;r2,S2) = J2 («^')n0m^np(^l)^n/i(^2)(-l)' 
nljm>0 

m(^25 ^ 2 ) 0 /jm(^ 2 : ^ 2 ) 0 /j—m(^l 7 ^l)] 

= E iuvU„^nnl,{n)nnlJ{r2){-ly-^+^Yl 

nljm>0 A,A' 

FA":i.yS;:i./<A(sJO«A.(sJ,) - c-y cy;i.yA(%)y,A.(ni) 
= E (™)..ii».’J»ij(’'i)K„y(''2)(-i)'“^+”E«A(ni)yiA'(n2) 

nljm>0 A,A' 

fxj—m _ fxjm ^j—m 

^lA,\silA',ls 2 ^lA',ls 2 ^lA,\si 


(A. 6 ) 


It is obvious that fi;(r,r', ),) 7 ^ —^(r, 4-; r', t); i-G- in the coordinate representation the spin 
structure of k has nothing common with (A.3). 
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The anomalous density defined by (A.6) has not definite angular momentum J and spin S. 
It can be represented as the sum of several terms with definite J, S. We have: 

0<J<2j 

^ ® 0i(2)}oo + X! ® 4>j{‘^)}j0- (A.7) 

We are interested in the monopole pairing only, so we omit all terms except the first one: 


® 0 ,( 2 )}oo 


2j + l 

Remembering the definition of 0 function we find 




(A.8) 


Sl)4>lj-m{^2, S2 )]j=q — ^ ^ 

A,aA',(T' 

The direct product of spin functions in this formula can be written as 


(A.9) 


Xi,(Sl)Xi^/(s2) = E'^fE.'{Xi(Sl) ® Xi(s2)}5S 


S',S 


— C'iE(7'{Ti('5i) ® yi (s2)}oO + E^|Eo-'{Ti('Sl) ® yi(s2)}lE- 


(A.IO) 


According to this result the formula for k consists of two terms: the one with S' = 0 and 
another one with S' = 1. It was shown in the paper [30] that the term with S' = 1 is an order 
of magnitude less than the term with S' = 0, so we can neglect by it. Then 


XiAsi)xia'is2) = (-1)5 (8)yi(s2)}oo 

1 1 

= (-l)5“'';^5a,-a' fyi i(si)yi 1 (ss) - X1 1 (si)y 1 1 (S 2 ) 


- 2 E-'t'(- 1 ) 


2 2 


S.iS 


1 


r) S‘2 ^ 5 l n S‘2 n 


(A.ll) 
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Inserting this result into (A.9) we find 


( 1 ) S2)]j^q — ^ Sg_^l5g^_l 


2 j + l 


A,A' l',cr 

~ 2 ~ 


2 j + l 


A,A' i^,tT Zt -\- I ’2 




21 + 1 


-ly-^ E >ZA(fIl)>;A'(fl2)(-l)X-A' 


A,A' 


^ k . -i,. .i„.l (-ir-'ipKcosni^), 


2 V -y (A.12) 

where PycosQu) is Legendre polynomial and f2i2 is the angle between vectors ri and r 2 . With 
the help of this result formula (A.6) is transformed into 

1 


K 


(ri,Si;r 2 ,S 2 ) 5 J = 




1 


51o^52-o Si-t^^S2 


2} 47]- 


E {uv)nljmT^nlj{ri)Tlnlj{r2)Pl{cOSni2^A.13) 


nljm>0 

Now it is obvious that in the coordinate representation k with J = 0,A = 0 has the spin 

structure similar to the one demonstrated by formula (A.3): 

0 K(ri,r 2 ) 

-fi;(ri,r 2 ) 0 


K 


(ri,Si;r2,S2)5J = 


(A.14) 


with 


1 

K(ri,r 2 ) = — E {uv)nijm'Pnij{ri)'R-nij{r2)Pi{cosQi2). 

^ nljm>0 


(A.15) 


Appendix B 


Wigner transformation 

The Wigner Transform (WT) of the single-particle operator matrix -Fri,o-;r 2 ,o-' is defined as 


[Fr 


r‘l5<T;r2,cr' 


WT = Fa 


.(r,p) = [ 


rH-s/2,cr;r—s/2,(j' 


(B.l) 


with r = (ri -|- r2)/2 and s = ri — r 2 . It is easy to derive a pair of useful relations. The first 
one is 


AA(r.P) = 


r+s/2,cr;r— s/2,(t' 


= / d'^se 


— 2 pS 


./np-- 


r—s/2,(j;rH-s/2,cr' 
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i.e., [-Fr\,„;r 2 ,<T']wT = [Frua'-,r 2 ,<T]wT = F*>^{r,p). The second relation is 

p) = F^<^/(r, -p) = J d^sF'^'^^^Fr+s/2,a-,r-s/2,a' 

For the hermitian operators p and h this latter relation gives [p*^ o-]wt = Pc 

fe,<7;r2,<7]wT = Ka{r,-p). 

The Wigner transform of the product of two matrices F and G is 

- - f ih 

[FG]wt = F{r, p) exp f — A j G(r, p), 

where the symbol A stands for the Poisson bracket operator 

-f— —> 

I (9 (9 d d 

^ 1 dvi dpi dpi dvi 


Appendix C 

Integrals of motion 

Isovector integrals of motion: 


2 8 


7^21 + 27^22 


+ V 3 + V 


1 


2 V^C2 

1 

\/3c2m 


(^(ci - C2)(ci + 2 c2) + 2ciC3 - 
^Ci + C 2 + 2 c 3 - VqI = const, 


tt 

00 




h 


Ct,-t'^FF 


— 3v^(l — q;)ko-R 2 o 


V3c2m v3c2 V o 


3 ~ Ao(r') J p 

ih-pc2Lu H-^— < in- 


^21 + ^(2ci + C2)7^2 i ~ Y -^^20 


~ (^^ 4 ^ ~ ~ ^ ^ 4:VSKoaLl^ nil 

— — 4c 2) — 2v^—Koct-^io^ "^ooj = const, 


(B.2) 

(B.3) 
(r, — p) and 

(B.4) 


const. 
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-pU 


22 


~ ("^20 + V^’^do) + y(Ci — C 2 ) VX^ — {T^X) + V^^do 

Ao(r' 


ih^iZ2i - H- 


^'O'l-t *^1 'O-l-t 

^ /Von r— /Von 


= const, 

2 pit 
' 00 


Q • “20 /7T ' “UU /7T 

O V'JC 2 y 6mC2 


■t,v^ Ao(r') 

‘'*2^^' - — 


/2^tt , 2 (ci + C2)^P ^ m (ci - C2)(ci + 2 c2)^P 
2' 20 + ^2^ ' 00 + ~ ^ '^1 


-^21 

where 


GhKoalCo 
Ao(r' 


' 2 ^tt *^1 Tptt 2 

Q '■^20 /77 '^00 /JT > 00 

O V'JC 2 y 6mC2 


h 


2 y3c2 

= const, 

'Pit 


00 




= const. 


= const. 


(C.l) 


Cl = 2ma;2 - Ca = 4:V6koRX^ + 

c^ = niu‘^ — 4A/3Q;fi;o-Roo + v^(l + a)KQR2o- 

Isoscalar integrals of motion are easily obtained from isovector ones by taking a = 1. In 
the case of harmonic oscillations all constants const are obviously equal to zero. 


Appendix D 

(D.l) 
(D.2) 

where x = 2app', 


n) = -Zp_p-<^p^ 




fi;'’(r,p') 0o(a^) ~ 4Q;V"^02(a^) e p’"^dp', 

J ^ -I 

J fi;’’(r,p')[0o(a^) — 2ap'^<pi{x)]e~°‘^'^p'^dp , 


1 1 
00 (a^) = — sinh(a;), 0i(a^) = 

Ob 


x^ 


cosh(a;)-sinh(a;) 


X 


x-^ 


3 \ 3 

1 H —-] sinh(a;)-cosh(a;) 


x^ 


X 


Anomalous density and semiclassical gap equation [17]: 

1 A(r,p) 


K(r,p) = 


/i 2 (r,p) + A 2 (r,p) 


(D.3) 


(D.4) 
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(D.5) 


/\ /|\ 

^'(Ip - P I) / 

y^/i2(r,p') + A2(r,p') 

where v(|p — p'|) = with (3 = —iVoKr^^Tr)^ and a = 
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